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A possibility to describe magnetism in the iron pnictide parent compounds in terms of the two- 
dimensional frustrated Heisenberg J1-J2 model has been actively discussed recently. However, recent 
neutron scattering data has shown that the pnictides have a relatively large spin wave dispersion in 
the direction perpendicular to the planes. This indicates that the third dimension is very important. 
Motivated by this observation we study the J 1 - J2- J c model that is the three dimensional generaliza- 
tion of the J1-J2 Heisenberg model for S = 1/2 and S = 1. Using self-consistent spin wave theory 
we present a detailed description of the staggered magnetization and magnetic excitations in the 
collinear state. We find that the introduction of the interlayer coupling J c suppresses the quantum 
fluctuations and strengthens the long range ordering. In the Ji-J2-J c model, we find two qualita- 
tively distinct scenarios for how the collinear phase becomes unstable upon increasing Ji. Eithcr 
the magnetization or one of the spin wave velocities vanishes. For S = 1/2 renormalization due to 
quantum fluctuations is significantly stronger than for S — 1, in particular close to the quantum 
phase transition. Our findings for the Ji-J2-J c model are of general theoretical interest, however, 
the results show that it is unlikely that the model is relevant to undoped pnictides. 



I. INTRODUCTION 

Over the last two decades, there has been considerable 
interest in the two dimensional (2D) S = 1/2 Heisen- 
berg antiferromagnet with frustrating interactions. One 
of the most widely studied models is the square lattice 
J1-J2 model, with both nearest neighbour J\ and second 
nearest neighbour J2 antiferromagnetic interactions. The 
Hamiltonian is 

h = Ji J2 Si • Sj + j 2 J2 s *- s j w 

<ij) ««» 

where stands for summation over nearest neigh- 

bors (NN) while stands for summation over next- 

nearest neighbors (NNN). Betts and Oitmaa^ were the 
first to point out that there is a finite long range order in 
the two dimensional Heisenberg model at zero tempera- 
ture. Early studies^ showed that the ground state of the 
pure Ji model has Neel order reduced by quantum fluctu- 
ations. The Neel order is destabilized with increasing J2 
and at some critical value of J2I J\ a phase transition to a 
quantum disordered phase occurs. On the other hand, for 
large J2 / J\ the system will order in a stripe-like fashion 
of alternating rows (or columns) of spins up and down. 
The long-range magnetic order is reduced by quantum 
fluctuations. As J%j J\ is reduced the collinear phase will 
become unstable at some critical ratio. There is substan- 
tial evidence; see Refs. HH! and references within, that 
the ground state of the quantum disordered phase has 
no long-range magnetic order and is dominated by short- 
range singlet (dimer) formation for 0.4 < J2/J1 < 0.6 for 
S = 1/2. The stability of such a configuration implies 
that the lattice symmetry is spontaneously broken and 
the ground state is fourfold degenerate. 



Using series expansion and mean-field spin wave the- 
ory methods, Singh et al. studied the excitation spectra 
of the square lattice J\- J2 Heisenberg antiferromagnet 6 . 
They showed the excitation spectra is gapless at only two 
symmetry related points of the Brillouin zone (0, 0) and 
(0,7r), whereas the accidental degeneracies at (7r,7r) and 
(ir, 0) are lifted by the 'order by disorder' effect, where 
the quantum fluctuations select a collinear ground state. 
Furthermore, they found the ratio of the spin-wave ve- 
locities along the x and y directions depends sharply on 
the J2/J1 ratio. 

Besides being of general theoretical interest the J\- 
J2 model is relevant to the real layered magnetic 
materials^. However, the real materials are not strictly 
two dimensional and contain a small interlayer coupling 
J c . For example, Rosner et. al. 9 found that J c /Ji ~ 0.07 
for Li 2 VOSi04, which can be described by a square lat- 
tice J1-J2 model with large J^^. This provides mo- 
tivation for studies of a three dimensional extension of 
the J1-J2 model. Such an extension for S = 1/2 has 
been recently studied using coupled-cluster and rotation- 
invariant Green's function methods^, a version of ef- 
fective field theoryii as well as different kinds of spin 
wave approache o 12 ' 13 . In particular these studies fo- 
cused on the influence of J c on the existence of an in- 
termediate quantum disordered phase at J2/J1 ~ 0.5. 
Schmalfufi et.alJ^ found that upon increasing the inter- 
layer coupling J c > the intermediate phase disappears 
at J c (0.4 — 0.6) J2. Our interest in the 3D model is 
mainly motivated by the discovery of superconductivity 
in the iron pnictides^. 

Parent pnictides demonstrate alternating spin stripes 
and therefore it is quite natural to assume that the 
collinear phase of J1-J2 model describes the system. 
Since their discovery, many investigations have been 
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FIG. 1. (color online) (a) Schematic diagram of the three di- 
mensional spin ordering in the Fe pnictides. Here we show the 
a — b plane for the considered model with nearest neighbour 
coupling Ji and next nearest neighbour coupling J2. The in- 
terplane coupling J c is directed into the page. In addition we 
show in (b) the real space positions of the lattice points A, a, 
B, and b which make up the unit cell. 



focused on understanding the magnetic properties of 
the pnictide parent compounds^—. Magnetic long 
range order has been established in LaOFeAs and 
Sr(Ba,Ca)Fe2As2 using neutron scattering^—, muon 
spin resonance (/iSR)^ and Mossbauer spectroscop y 29 ' 30 . 
The neutron studies reveal the parent compounds dis- 
play a columnar antiferromagnetic ordering with a stag- 
gered magnetic moment of (0.3 — 0.4)//^ in LaOFeAs and 
(0.8 — 1.01)/kb in Sr(Ba,Ca)Fe2As2. In this columnar ar- 
rangement, stripes of parallel spin order along the b axis 
and antiferromagnetically along the a and c axe a 27 ' 31 ; see 
Fig.ffl 

The spin wave velocities along the a and c axes 
have been measured in neutron scattering studies of 
SrFe 2 As 2 2 3 - and BaFe 2 As2^ with v a « 205 meV and v c s» 
45 meV. Upon lowering the temperature the parent pnic- 
tide compounds undergo a structural lattice distortion 
from a tetragonal to an orthorhombic structure. The or- 
thorhombic distortion is very small, and is of order a frac- 
tion of one per cent 3 - 2 -. The structural transition happens 
at a temperature slightly higher or equal to the mag- 
netic ordering temperature. This coincidence naturally 
suggests that the structural transition can be driven by 
the spontaneous violation of the Z2 symmetry (nematic 
transition) in the spin stripe phas o 15 ' 16 . Note that this 
suggestion is based purely on symmetry arguments and is 
valid for both the "localized" and "itinerant" paradigm. 
Let us call this scenario, scenario A. The structural tran- 



sition can also be fully independent of the spin structure 
(scenario B). Moreover, the transition can drive the spin 
structure (scenario C). In the present work we ignore 
the very small orthorhombic distortion and consider the 
"isotropic" model i.e. a tetragonal lattice. In relation to 
the iron pnictides this approach makes sense in cases A 
and B, and it is not justified in case C. 

In the tetragonal phase, the Fe sites form square pla- 
nar arrays, such that the sites of adjacent planes lie above 
each other; see Fig. [T] Recently, there has been debate 
concerning the spin wave velocity along the 6-axis i.e. 
along the spin stripe s 23 ' 33 . On the one hand direct neu- 
tron scattering data from twinned samples 23 indicate the 
value of Vb is comparable with v a . On the other hand 
analysis^ 3 - of the NMR relaxation rate indicates an order 
of magnitude smaller value of w&, Vb ~ 10-30 meV. The 
small value of Vb implies that the system is close to a 
quantum critical point. Therefore from now on we will 
refer to the small Vb scenario as the critical one and the 
large Vb scenario as the non-critical scanario. 

Band structure calculation s 34 ' 35 have shown that 
columnar antiferromagnetic ordering is the most stable 
structure. On the one hand, band structure results in- 
dicate a local moment of up to 2.3/ie per Fe sit o 17 ' 35 ' 36 . 
This value is too large compared to experiment. This 
has led to the suggestion that the ordered moment might 
be strongly renormalized by magnetic fluctuations, as de- 
scribed by the frustrated 2D Ji- J 2 modeli&Si 3 ^ 3 - 7 -. This 
suggestion implies the critical scenario because of very 
strong quantum fluctuations. On the other hand, there 
are studies suggesting that the smallness of the staggered 
magnetic moment can be explained by electronic effects 
such as hybridization and spin-orbit coupling 22,38 . 

Motivated by what is currently known about the pnic- 
tides the goal of the present work is three- fold. First, 
we provide a quantitative theory for the magnetic exci- 
tations based on a minimal spin model, namely the J\-J% 
Heisenberg model with interlayer coupling J c for S — 1/2 
and S — 1. In particular, we will show the interlayer 
exchange coupling suppresses the quantum fluctuations, 
strengthens the staggered moment, and hence dramati- 
cally increases the stability of the columnar phase. Fur- 
thermore, it has been recently suggested that the strong 
reduction of the magnetic moment as seen in 2D is not 
possible in 3D 3 - 9 - because the 3D coupling cuts the loga- 
rithmic divergence of quantum fiucuations. In this paper 
we will show that a considerable reduction of the stag- 
gered moment is still possible for S — 1/2 for small values 
of J c while a significant renormalization for S = 1 is only 
possible for extreme fine-tuning. 

Second, we discuss the phase diagram and the exis- 
tence of a dimerized phase for the J1-J2- J c Heisenberg 
model. Indeed we find a dependence of the point where 
the staggered magnetization vanishes as a function of in- 
terlayer coupling J c . In particular, two different insta- 
bilities of the columnar phase appear: Either the long- 
range order or one of the spin wave velocities vanishes. 
Note that in three dimensions the vanishing of one of 
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the spin wave velocities does not imply the vanishing of 
the staggered magnetization. We will discuss these two 
kinds of instabilities in a forthcoming publication in more 
detail^. 

Third, we make quantitative predictions of the disper- 
sion and spin wave velocities parallel (v a ), perpendicu- 
lar (vb), and through (i; c ) the magnetic stripes; the ratio 
Vb /v a is a useful way to determine the degree of magnetic 
frustration Jx/J%- The spin wave spectra over the full 
Brillouin zone show dramatic differences between a sys- 
tem deep in the columnar phase and one close to a quan- 
tum critical point. This provides a robust experimental 
way to distinguish the critical and non-critical scenarios 
presented above. In addition, we will compare our model 
with experimental evidence and discuss whether the J\- 
Ji-3 C is an appropriate model for describing magnetism 
in the pnictide parent compounds. 

II. MODEL AND METHODS 

It is well known that the parent iron pnictidcs arc not 
simple Mott-Hubbard insulators, but rather are bad met- 
als with a very small Fermi surface. On the one hand this 
supports an itinerant picture for the compounds while on 
the other hand there are many experimental indications 
for localized moments. On a pure theoretical level it is 
clear that the Hubbard on-site repulsion is rather large 
and therefore even an itinerant system must be close to 
the Mott-Hubbard regime. Most likely the truth is some- 
where in the middle since there are localized and delo- 
calized degrees of freedom, and it is not clear yet how to 
combine these two descriptions. 

In the present work we disregard the itinerant degrees 
of freedom and consider a model of well-localized spins 
described by a Heisenberg Hamiltonian consisting of ef- 
fective in-plane nearest neighbour J\ , next nearest neigh- 
bour Ji and interlayer J c exchange interactions 

H = J\ Si ■ Sj + J2 Si ■ Sj + J c Si ■ Sj (2) 
where and correspond to summation over 



nearest-neighbour (NN) and next-nearest neighbour 
(NNN) pairs in the plane and \i,j \ corresponds to nearest 
neighbour pairs between the planes, see Fig. Q] 

Since we want to describe the observed spin-stripe 
phase in the parent pnictides we are interested in the re- 
gion of parameters of the effective Heisenberg model that 
supports the phase. It is well known that this phase is 
stable at sufficiently large J2, J2 > 0.5Ji. The particular 
values of J\ , J2 , J c can be determined from comparison 
with known magnetic excitation spectra and we will dis- 
cuss this issue later in the paper. Note that we consider 
the "isotropic" Heisenberg model, the value of J\ is the 
same in both directions, along the spin stripes and per- 
pendicular to the stripes so that we consider the entire 
stripe ordering to be purely spontaneous. 

We calculate the sublattice magnetization and mag- 
netic excitations using self-consistent spin wave theory 
for 5 — 1/2 and 5 = 1. This method has already been 
shown to work well for 5 = 1/2 and 5 = 1 in the colum- 
nar phase for the two dimensional cas o 6 ' 18 i 41 . We have 
used the Dyson-Malee v 42 i 43 as well as the Schwinger bo- 
son representation— which on the level of self-consistent 
mean field theory yield the same results for T — 0. The 
boson operators Ai, a^, Bj, bj on the respective sub- 
lattices, see Fig. [TJ are introduced in the usual way by 
performing a Dyson-Maleev^ 2 - transformation of the spin 
operators. 

Sf = S - a^ai 

Sj = (25 - ajai)a,i 

S i = a i f 

S* = -S + b/b 3 

5t = -6/(25-6A) 
57 = -6, (3) 

The definition of A and B is similar to a and b. Using Q 
the Hamiltonian ([2]) may be presented to quartic order 
in the operators A, a, B and b 
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H = H Aa + H Bb + H Ab + H aB + H ab + H AB (4) 



H Bb = 



H 



Ab 



H 



uB 



ab 



-JiS 2 + Ji ^ ls(A l ^A l + a/ajj - Aidj - A^a/) + ^(A^AiAidj - 2A i i A i a j ' i aj + A^a^a^aj) 1 

- J^ 2 + JiJ2 \s(BjBi + bjhj - B t b, - B s: V) + \{BjBiBibj - 211/ II, h/!^ + Bjbj%%)\ 

JiS 2 -JiJ2 \S(AMi + a^aj - A^ - Ajbj) + \{A/ A^b, - 2AM i b j h j + A/h/b/hA 

JiS 2 -Ji^2\ S(aMi + B^ Bj - diBji - ajBj) + hd/a^Bj - 2d { •,,!!/ IS, + a^B^B^Bj) 1 
(id) 1 J 

H ab = -J 2 S 2 + J 2 ^2 | S(a^ai + b^bj - a,ibj - ajb^) + -(a, 1 a^bj - 2a l t a i 6/6 J - + a^b^b^bj) i 
({id)) ^ ' 

- J C S 2 + Jc^2) S(a^di + b^b-j - aibj - ajbj*) + ^-(d^didibj - 2u/ ti,h/h, + a/bj^b^bj) 1 
Hd\ 1 J 

// lB = -J 2 S 2 + J 2 \s(AMi + B^Bj - AiB-j - A^B/) + 1(A/ A^Bj - 2AM i B^B j + A, 11/ 11/ II, A 
({id)) L j 

J C S 2 + J c ls(A/A l + Bj^Bj - AiBj - AjB*) + ^(AMiAiBj - 2A/ A, IS/IS, + A/ 11/ II, II, A . 

\ij\ L " J 



We perform the Hartree-Fock mean field decoupling of formed to 
quartic terms using the following notations. 

/ - ~ = (AM,) = (of*) = (Bj^Bj) = (6/6,.) 

F = (Aibj) = (Ai%*) = (enBj) = (a/B^) 
G = (A/ aj ) = (Aid}) = (Bi\) = (Bib}) 

g = (A t Bj) = (AjB*) = (aibj) = (a, V) 

h = (AiBj) = (AjBS) = (dibj) = (djbl) . (5) 

The difference between q and h in (TST) is that q corre- 

AF1 7"F AF2 

sponds to the expectation value of next nearest neighbour H = a + H + H + H 

pairs in the plane and h corresponds to nearest neighbour a — —2J 2 S 2 — 2J C S 2 

pairs between the planes, see Fig. [1] It is convenient to „ afi _ 7 c v~V a t a t into l u 
. , jjcj li — Ji/io J (Ai Ai + a,j dj + r>i r>i + Oj Oj 

introduce parameters v, rj and % denned as ^ — ' 

(i<j) 

- Aidj - A/a^ - Bibj - Bi t 6 J - t ) 
H F = -JtvS^iAM, + bjhj + djdi + II II 

(i,3) 

- AJ,, - Ai% - o.li: - a, IS,, 

H AF2 = -J 2 rjS ^2 (a^di + bjhj + AMt + B^Bj 

((id)) 

- dibj - ajbji - AiBj - AjBj) 

Values of the parameters obtained in the self-consistent + JcXS^2i a i a t + bj bj + A % A % + Bj Bj d % bj 

procedure described below are plotted in Fig. [5] After the ' lJ 

Hartree-Fock decoupling the Hamiltonian ^ is trans- — ajbj' — AiBj — A^Bj^) (7) 



fiS = 


S^ 


1 

" 2 


-/- 


-G 


vS = 


s^ 


1 

" 2 


-/- 


-F 


V S = 


s^ 


1 

" 2 


-/- 


"5 


\S = 


s^ 


1 

" 2 


-/- 





(0) 
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In the momentum representation the Hamiltonian reads 



H = 



H 



a 

AFl 



H 



AF2 



H AF1 + H AF2 
2 



a + H 

-2(2J 2 + J C )NS 



flfctflfc + B^B k + b k %) 

fcGMBZ 

C-(A k a- k + i4* t o_ fc t + B fc 6-fc + B^b^)) 
-AJivS ((A fc t A fc + a fe t a A; +B fc t B fe + 6 fc t 6 fe ) 



£T = 



E 

fcGMBZ 



C*+(A fe & fe t + Afetfofe + (a k B k f + a^Bk)) 
45(2 J 2 r) + J cX ) V ((AfcUfc + a fc ta fc 



E 

fcGMBZ 



+ B k ^B k + b k f b k ) - fik(A k B- k + AjB.J 



(8) 



Because of four different sublattices the Fourier trans- 
form is defined in Magnetic Brillouin Zone (MBZ), this 
is the Brillouin Zone of one sublattice, say the sublatticc 

a, 



1.5 



0.5 





S = 1/2^—- 











0.5 



1.5 



1.5 




E 

fcGMBZ 



(9) 



In Eq. ((5]) H F and H AF1 consists of ferromagnetic and 
antiferromagnetic intraplane terms for nearest neighbour 
pairs (Ji) while H AF2 consists of antiferromagnetic in- 
traplane next nearest neighbour pairs (J 2 ) and antiferro- 
magnetic interplane nearest neighbour pairs ( J c ). In ad- 
dition for simplicity we have introduced the coefficients 
C + , C_ and [if. 



J 2 r}(cos(k x ) + cos(fc y )) + J c x cos(fc z ) 

2 J 2 r) + J c x 
cos((k x + k y )/2) = cos(fch/v / 2) 
cos((k x - ky/2) 



Mfe 

C+ 

C 2 = cos(fc z ) = cos(fc c ) 



COS 



(k a /V2) 



(10) 



The components k x and k y are directed along diagonals of 
the base square, see, Fig. [1] At this stage it is convenient 
to unfold the Magnetic Brillouin Zone and to use the full 
Brillouin Zone. 

k = (k ai fc;,, k c ) 

— 7T <k a < 7T 

— 7T <kb< 7T 

-7r<A: c <7r (11) 

Therefore the momentum summation in all subsequent 
equations is defined as 



E 



N 
~2 



dk a 
2tt 



dkj, 



dk c 



(12) 



Since we are considering the four distinct sublattices A, 
a, B, and b the full Brillouin zone over counts the num- 
ber of degrees of freedom. To compensate this we have 



1.5 



0.5 



1.5 



0.5 



0.5 



1 

=vj 2 



0.5 



0.5 



1 



1.5 




1.5 




1 

■hi*! 



1.5 



FIG. 2. (color online) Quantum correction parameters /i, v, 
f/, X as a function of J1/J2. Calculations were performed by 
numerical iteration for a range of y values for both 5 = 1/2 
and S = 1. Here we show the results for y = 0.01 (solid) and 
0.10 (dashed). 
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introduced an additional prefactor in the definition of 
summation (|12l) . 

In the symmetry broken phase the dispersion reads 



(13) 



The coefficients Ak and By. are easily obtained by di- 
agonalizing the Hamiltonian (J8| in the usual way via a 
Bogoliubov transformation 

A k = 2J 2 S{\ + xvC+) 

B k = 2J 2 S{2 V C+C- + xiaC- + y X C z ) (14) 

where for simplicity we have introduced the frustration 
parameters 



x = Ji/J-2 

y = J c /J 2 
and have the numerical factor 

\ = x(fi - v) + 2rj + yx 



(15) 



(16) 



which is expressed in terms of the quantum correction pa- 
rameters ([B]) . The sublattice magnetisation m s is defined 

as m s — (Si Z ) so that 



N 4^ 



.4, 



VA k 2 - B k 2 



(17) 



The quantum correction parameters were determined 
from the self consistent equations by numerical iteration; 
see Eq. (0 and Eq. (J6j) . 



F 



N ^ 
k 



A k 



Al B k 2 
BkC- 



G 

9 
h 



/. V Ap. 2 — By, 2 
1 AkC + 
N ^ 



k yAk z — Bk 



n 4^ 



2 _ p j 

Bkf^k 



k yAk—Bk 
1 ^ BkC z 
N ^ 



k V Ak 2 — Bk 



(18) 



For 2 J 2 > J\ 3> J c > 0, the classical ground state is 
the stripe ordered phase with ordering vector (7r,0,7r). 
In addition to breaking spin rotational and time-reversal 
symmetries, the lattice symmetry is spontaneously bro- 
ken, the ordering vector (0, tt, tt) is equally possible. To 
be specific we consider the (tt, 0, it) phase. The excitation 
spectrum is gapless at (0,0,0) and (tt, 0,tt). These are 
the Goldstone modes with dispersion 
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S= 1/2 
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FIG. 3. (colour online) Staggered magnetization in the three- 
dimensional magnetic Brillouin zone as a function of the ratio 
of the exchange couplings Ji and J2 . The calculated three di- 
mensional staggered magnetisations were obtained using (a) 
linear spin wave theory and (b) self consistent spin wave the- 
ory. In each plot we show results for 5 — 1/2 and 5 = 1 for 
y = 0.01 (solid) and 0.10 (dashed). 
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Columnar spin-stripe , 



Columnar spin-dimer 



w(q) ~ V Va 2 q a 2 + v b 2 q b 2 + v c 2 q c 



(19) 



1.8 1.85 1.9 1.95 

X = V J 2 



FIG. 4. (colour online) The phase diagram of the 5 = 1/2 
3D J1-J2-JC Heisenberg model, y = J c / J2 measures the inter- 
layer coupling. The region above the curve corresponds to the 
columnar spin-stripe phase and the region below to possibly 
the columnar spin-dimerized phase. 
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FIG. 5. (colour online) The phase diagram of the 5 = 1 3D 
Ji-Ji-Jc Heisenberg model, y = J c jJi measures the inter- 
layer coupling. The region above the curve corresponds to 
the columnar spin-stripe phase and the region below to pos- 
sibly either a columnar spin-dimerized phase like for 5=1/2 
or a quadrupoler ordered phase. Note the very different scales 
compared to Fig. [4] 



where q is measured relative to the ordering vector and 
v a , Vb, v c are the spin wave velocities along the crystal 
axes 



v a = 2J 2 S y/ (277 + EM) ( 2r 7 + xfi + yx) 
Vb = 2J 2 Sy / (2i] ~ xv)(2r\ + xfi + yx) 



r, 



= 2 J 2 S y/ yxi^V + + yx) 



(20) 



The existence of this branch of magnetic excitations fol- 
lows from the Goldstone theorem. There exists also a 
second branch of spin waves for cj(0,7t,0) = w(7t,7t, w) 
and w(0,7r, ir) — w(7r,7r,0) which have not yet been ob- 
served experimentally. Since they are gapped the inten- 
sity of scattered neutrons would be substantially lower 
than the primary branch at (jr, 0,tt). Their energies read 



w(Q,7r,0) 

w(0,7T,7t) 



4J 2 S\/(2r] - xv){x[i - xv + yx) 
4J 2 Sy/(xn - xv)(2r] - xv + y X ) (21) 




> 





(c) 









^^^^ S^l\^- 

x 




s = 1/2V \ 
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X = Jj/Jj 



III. GENERAL RESULTS 

The influence of the ratio of exchange parameters 
y = J c / ' J 2 on the resulting staggered magnetization, spin 
wave velocities and excitation spectra is studied for a 
range of values of y. The staggered magnetisation was 
calculated by numerical iteration for a range of y assum- 
ing y <C 1; the results are shown in Fig. [3l We show re- 
sults for both linear spin wave theory and self -consistent 
spin wave theory for S = 1/2 and 5=1. 

We have also calculated the "critical" value of x — 
J\j J 2 where either the staggered magnetization or one 
of the spin wave velocities vanishes, both, for S = 1/2 
in Fig. @] and for S = 1 in Fig. \5\ A detailed distinction 



FIG. 6. (colour online) Spin wave velocities along the crystal 
axes as functions of the ratio J1/J2 for y — 0.01 (solid) and 
0.10 (dashed). Results are presented for 5 = 1/2 and 5=1. 
Note that the ordinate values are divided by y/y for v c only. 

between vanishing staggered magnetization or vanishing 
spin wave velocity will be presented elsewhere^. Basi- 
cally these figures show the phase diagram of the model 
as it results from self-consistent spin wave theory. For 
S = 1/2, the regions above the lines correspond to the 
columnar spin-stripe phase and the regions below them 
correspond to the columnar spin-dimerized phase. It is 
well known that the transition between these phases is 
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of first order. Therefore, strictly speaking, the criterion 
of a vanishing magnetization is not quite the correct one 
to indicate the transition, the true transition happens at 
a slightly smaller value of J1/J2 than that indicated in 
Fig. |4j However, it is known, see e.g. Ref. H, that the 
criterion gives practically the correct value of the critical 
point for y = 0. Here we assume that the same is true 
for small nonzero y. 

In Fig. [3] we observe an interesting behaviour of m a (x) 
as a function of increasing interlayer coupling. In the 
vicinity of x — 2 and for small values of y the renor- 
malized magnetization changes very rapidly on small pa- 
rameter changes, differing significantly from linear spin 
wave theory calculations even for 5=1. This clearly in- 
dicates strong quantum fluctuations. These fluctuations 
are stronger for 5 = 1/2, see for instance the larger de- 
viation of the region of instability from the value x — 2, 
than for 5 = 1 (compare Fig. @] and Fig. [3]). This also 
implies that very small values of the staggered magneti- 
zation can be obtained for 5 = 1/2 more easily, i.e., with 
less fine-tuning, than for 5 = 1. But qualitatively, the 
curves for 5 = 1/2 and for 5 = 1 are very similar. For 
all parameters, the long-range order is strengthened by 
the coupling in the third dimension. 

We would like to briefly comment on the intermedi- 
ate quantum phases for 5 = 1/2 and 5 = 1 near x = 2. 
There is a consensus that in the 5=1/2 two-dimensional 
J1-J2 model there is an intermediate magnetically disor- 
dered phase at OAJi < J 2 < O.6J1. We believe that this 
is the columnar spin dimer-phase. According to our data 
for the 3D model presented in Fig. 0] the columnar spin 
dimer phase disappears around J c j Ji ~ 0.25. The criti- 
cal value of J c for disappearance of the magnetically dis- 
ordered phase differs from that determined previously-^ 
by the couple-cluster method, J2/J1 ~ 0.36. We do not 
think that the difference in the value is significant, both 
our method and the method of Ref. [l^ are approximate. 
More importantly, there is qualitative agreement about 
the phase diagram. For 5 = 1 we also found a tiny 
region of an intermediate nonmagnetic phase shown in 
Fig. [5] This is qualitatively different from previous stud- 
ies^ that were unable to identify an intermediate phase. 
Note the very small scale on which we find the instability 
of the columnar stripe order. Unfortunately within the 
present method we cannot determine the exact nature of 
the phase and therefore in Fig. [S] the phase is shown by 
a question mark. 

The ratio of the spin wave velocities along the differ- 
ent directions is very sensitive to the value x = J\j Ji- 
The point is that in a real compound the staggered mag- 
netization can depend on a range of additional uncon- 
trolled variables such as itinerancy, hybridization, etc. 
The dispersion relation however only depends on the ef- 
fective Hamiltonian and thus is less ambiguous^. Our 
results for the spin wave velocities along the three crys- 
tal axes are shown in Fig. [5] versus the ratio J1/J2 for 
y = 0.01, and 0.10. The velocities are given in units of 
J 2 ; v a and only weakly depend on y, while on the 
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FIG. 7. (colour online) Spin excitation spectra along high- 
symmetry cuts through the Brillouin zone for 5=1. We 
show the excitation spectra for both a system deep in the 
columnar phase (J1/J2 = 1) and one near the quantum phase 
transition (J1/J2 = 1.98). Calculations were performed using 
self-consistent spin wave theory and here we compare results 
for y — 0.01 (solid and dotted) and 0.10 (dashed and dot- 
dashed) . 



other hand v c oc ^fy. The dependence of the ratios of the 
spin wave velocities on the values of the exchange cou- 
plings is stronger than the corresponding dependence of 
the staggered magnetization. So it is more appropriate 
to determine the values of the couplings from the spin 
wave velocities. 

Spin excitation spectra for 5 = 1 along high-symmetry 
cuts through the Brillouin zone are shown in Fig. [7l We 
show the excitation spectra for both a system deep in the 
columnar phase (J1/J2 = 1) and one near the quantum 
phase transition (J1/J2 = 1.98). In the present work we 
analyze the dependence of the spin wave dispersion on the 
interlayer coupling J c . The case of the small coupling, 
y = J C J J2 -C 1 is of special interest. Expanding the 
spectrum as defined in (|13[) in powers of y we find 

w 2 (k) =wj}(k) +y Suj 2 (k) (22) 
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Material 


LUA Moment (^b) 


Expt. Moment (/is) 


j , . ( \ T7U A o 

YjQXj r e/YS 




U.oO— 


NdOFeAs 


1.49 


0.25^- 


CaFe 2 As 2 


1.51 


0.80^ 


BaFe 2 As 2 


1.68 


0.87^. 


SrFe 2 As 2 


1.69 


1.01^ 



TABLE I. Comparison of the magnetic moment in units of 
/is as predicted from LDA calculations and those observed in 
experiment. 

where 

Lu 2 (k) = 4,JiS 2 [(x(fi -v)+2r] + xvC+f 

- (2r]C + C- +x[iC-) 2 } (23) 

and 

6io 2 (k) = &J%S 2 x[(x(ii ~v) + 2r,+ xvC+) 

- (27/C+C- + xfiC-)C z ] (24) 

Note that all parameters in Eq. (|2"3")1 and Eq. (|2"4")l are 
calculated at y = 0. Eq. (1231) is the dispersion in the 
2D case, it has Goldstone modes (k a ,kb) = (0,0) and 
(k a ,kb) — (tt, 0). The expanded 3D dispersion (|2"2")l 
has the Goldstone modes at (k a ,kb,k c ) = (0,0,0) and 
(k a ,kb,k c ) — (tt, 0, 7r) as expected. Note that w(k) is 
generally a non-analytic function of y at small y while 
w 2 (k) is the analytic one. This is why the expansion 
(f2"2"|) is written in terms of lu 2 . In the limit y — > the 
parameters [i,v,r\,x depend on i — J1/J2 only. The 
corresponding plots are presented in Fig. [2J These plots 
together with Eq. (|2"21 allow one to determine the spin- 
wave dispersion at arbitrary small y. 

IV. APPLICATION TO IRON PNICTIDES 

The smallness of the measured magnetic moment in 
the iron pnictides relative to theoretical calculations is a 
matter of controversy; there exist two different scenar- 
ios which offer different explanations for the discrepancy 
between experiment and theory. It has been suggested 
that magnetic fluctuations may strongly reduce the local 
magnetic moment, with the ratio x = J1/J2 being fixed 
to an appropriate value in the critical scenario. The al- 
ternative derives its explanation from the role of the local 
electronic orbitals and therefore the magnetic couplings 
are not determined by the value of the magnetic moment. 
Band structure calculations have shown that J\ and J2 
are antiferromagnetic and very similar in valu o 17 ' 35 i 36 . In 
Table U we compare the size of the experimentally mea- 
sured magnetic moment to that calculated by LDA meth- 
ods. Typically the experimental values are at least twice 
smaller than the LDA values. Due to strong quantum 
fluctuations in the model in principle one can obtain the 
required suppression of the staggered magnetization by 
two times by choosing x « 1.99 for S = 1, see Fig. [3Jb). 




FIG. 8. (colour online) Spin wave dispersion in the plane 
of the spin stripes. The spin wave dispersions u;(k) at 
(jj(k a ,kb,n) in units of J 2 is shown for Ji/J 2 = 0.76 (top) 
and J1/J2 = 1.972 (bottom) for y — 0.10. Values were ob- 
tained using self consistent spin wave theory for S = 1. 



However, due to the large change in m s upon small pa- 
rameter changes, it is clear from Fig. [2[b) that consider- 
able fine tuning is required. 

A more effective probe for the value of the ratio x = 
J\ j J2 is considering the spin wave velocities since they 
depend on the Hamiltonian only. In the present work we 
compare the J1-J2- J c Heisenberg model to experimental 
evidence; this is studied via the critical and non-critical 
scenarios. The dispersion in the plane of the spin stripes 
uj(k a , fcb, tt) is shown in Fig. [5J For the critical scenario 
we adjusted the ratio x = J1/J2 = 1.974 to fit the ra- 
tios Vb/v a and v c /v a obtained from the analysis^ 3 - of the 
NMR relaxation rate. For the non-critical scenario we 
used the experimental ratios of the spin wave velocities 
Vb/v a ~ 0.70 and v c /v a « 0.25 that follow from inelastic 
neutron scattering data^. In this case x — 0.76 gives 
the best fit. In addition we plot in Fig. [9] the dispersion 
of spin waves in the plane along and through the spin 
stripes ui(tt, k^, k c ). The dispersion uj(tt, fcf,, k c ) is helpful 
since it shows the effect of the interlayer coupling and 
how extensive the spin waves propagate along the c-axes. 
In addition it is a way to determine if strong quantum 
fluctuations persist in three dimensions. 

The plots in Fig. [S] and Fig. [S] provide a general quali- 
tative overview. A quantitative comparison between the 
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M{7T,kb,k c )/J 2 2 




FIG. 9. (colour online) Spin wave dispersion in the plane 
along and through the spin stripes. The spin wave dispersions 
oj(k) at uj(n,kb,k c ) in units of J2 is shown for J1/J2 = 0.76 
(top) and J1/J2 = 1.972 (bottom) for y = 0.10. Values were 
obtained using self consistent spin wave theory for S = 1. 



critical and noncritical scenarios and the inelastic neu- 
tron scattering data is depicted in Fig. [TO] We compare 
the fitted dispersions for the critical and noncritical sce- 
narios. The noncritical scenario agrees nicely with the 
experimental data whereas the critcal scenario does not 
describe the dispersion at higher energies due to problems 
with the reduced staggered moment and the anisotropic 
spin wave velocities. Therefore, it is quite clear that the 
pnictides, if decribed by the J\-Ji-J c Heisenberg model, 
are deep in the columnar phase with Ji = O.76J2. 

In addition to determining if the pnictides are in the 
critical or non-critical regime we study whether or not 
the J\-J2-J c Heisenberg model is an appropriate model 
for describing the undoped pnictides. In the critical sce- 
nario, Vb is small compared to v a implying a high density 
of magnetic excitations^ while in the non-critical sce- 
nario Vb is sizeable compared to v a implying a low den- 
sity of magnetic excitations. In Fig. [10] we observe that 
while the noncritical scenario does reproduce the known 
spin wave velocities, the dispersion cannot be matched 
globally. Significant differences persist at intermediate 
energies for the (n,kb,ir) dispersion curve. Therefore, it 
is clear that the J\-J2-J c Heisenberg model is not con- 
sistent with the data by Zhao et al^. We stress that 
the data is not consistent with the "isotropic" Heisen- 
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FIG. 10. (colour online) Comparison of fitted dispersions 
for the critical and noncritical scenarios with inelastic neu- 
tron scattering data. We see that the noncritical scenario 
agrees nicely with the experimental data whereas the critcal 
scenario does not describe the disperison at higher energies 
due to problems with the reduced staggered moment and the 
anisotropic spin wave velocities. 



berg model considered in the present work because of the 
minute orthorhombic distortion. On the other hand, the 
data is consistent with the anisotropic Heisenberg model, 
see discussion in Ref. online^. The anisotropic Heisen- 
berg model implies that there are additional degrees of 
freedom (orbital?) and it is outside of the scope of the 
present work. 

It has been recently suggested that the strong re- 
duction of the magnetic moment possible in two 
dimensions^ is not possible for substantial three di- 
mensional coupling and/or magnetic anisotropy because 
these additional couplings dramatically suppress the 
quantum corrections to the ordered moment^. The log- 
arithmic divergence of the quantum corrections to the 
staggered magnetisation seen in the square lattice J1-J2 
model for x — > x c is cut off by the addition of either an 
anisotropy gap or a third dimension. In Fig. [3] we do in- 
deed see an enhancement of the staggered magnetisation 
compared to the two dimensional case^ for increasing 
values of y for 5 = 1 /2 and 5 = 1. 

Since the quantum fluctuations are more significant for 
5 = 1/2 than for 5 = 1 their suppression by the three di- 
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mensional coupling is seen more distinctly in the 5=1/2 
data than in the 5=1 data. In particular, close to 
the critical values x c , where the columnar striped phase 
becomes unstable due to quantum fluctuations, the sup- 
pression for increasing three-dimensional coupling y — > 1 
is most clearly seen, cf. Fig. [3j which is in accordance 
with the findings by Smerald and Shannon 3 -. The quan- 
titative difference between the results for 5=1/2 and for 
5 = 1 is most strikingly seen in the difference of scales 
of the Figs. H and EJ 

V. CONCLUSIONS 

We have provided a quantitative theory for the mag- 
netic excitations on a tetragonal lattice based on a 
minimal spin model, namely the J\-J<2,-J c Heisenberg 
model. First, we have shown that the columnar phase 
is stabilized by the introduction of the interlayer cou- 
pling J c . Such a three-dimensional coupling dramatically 
strengthens the staggered magnetization and suppresses 
the strong quantum fluctuations. Since for the 5 = 1/2 
case the staggered magnetization is more strongly renor- 



malized by the quantum fluctuations than for 5 = 1 the 
effects of the suppression of the quantum fluctuations are 
more clearly seen for 5=1/2 than for 5 = 1. Yet both 
spin species behave qualitatively the same 

In addition, we have also shown that the position of 
the critical point depends on the value of of the relative 
interlayer coupling y. Again, this dependence is more sig- 
nificant for 5=1/2 than for 5=1 because the influence 
of quantum fluctuations and thus also of their suppres- 
sion is stronger for smaller spin than for larger spin. As 
function of y we found to distinct ways how the colum- 
nar phase becomes unstable. Either the magnetization 
or one of the spin wave velocities vanishes. 

Second, the strong reduction of the magnetic moment 
is possible in three dimensions if one considers small val- 
ues of y. One has to approach x c very closely even for 
5 = 1/2. Third, we have shown that one can conve- 
niently model uJk(y) for small y analytically. 

Finally, comparing the J\-Ji-Jc Heisenberg model with 
experimental data we found that such a model does not 
explain the data3^. If reproduces the spin wave velocities, 
it cannot match the dispersion globally. It is clear that 
further work is called for. 
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